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Background: The ground state neutron density of a medium mass nucleus contains fundamental 
nuclear structure information and is at present relatively poorly known. 

Purpose: We explore if parity violating elastic electron scattering can provide a feasible and model 
independent way to determine not just the neutron radius but the full radial shape of the neutron 
density Pn{r) and the weak charge density pw{r) of a nucleus. 

Methods: We expand the weak charge density of "^^Ca in a model independent Fourier Bessel series 
and calculate the statistical errors in the individual coefficients that might be obtainable in a model 
parity violating electron scattering experiment. 

Results: We find that it is feasible to determine roughly six Fourier Bessel coefficients of the weak 
charge density of '^^Ca within a reasonable amount of beam time. However, it would likely be much 
harder to determine the full weak density of a significantly heavier nucleus such as ^®®Pb. 
Conclusions: Parity violating elastic electron scattering can determine the full weak charge density 
of a medium mass nucleus in a model independent way. This weak density contains fundamental 
information on the size, surface thickness, shell oscillations, and saturation density of the neutron 
distribution in a nucleus. The measured pwir), combined with the previously known charge density 
pch(r), will literally provide a detailed textbook picture of where the neutrons and protons are 
located in an atomic nucleus. 

PACS numbers: 25.30.Bf, 27.40.+z, 21.10.Gv, 21.10.Ft 


I. INTRODUCTION 

Where are the protons located in an atomic nucleus? 
Historically, charge densities from elastic electron scat¬ 
tering have provided accurate and model independent 
information [T]. These densities are, quite literally, our 
picture of the nucleus and have had an enormous impact. 
They have helped reveal the size, surface thickness, shell 
structure, and saturation density of nuclei. 

Where are the neutrons located in an atomic nucleus? 
Additional, very fundamental, nuclear structure informa¬ 
tion could be extracted if we also had accurate neutron 
densities. For example, knowing both the proton and the 
neutron densities would provide constraints on the isovec¬ 
tor channel of the nuclear effective interaction, which is 
essential for the structure of very neutron rich exotic nu¬ 
clei. 

However, compared to charge densities, our present 
knowledge of neutron densities is relatively poor and may 
be model dependent. Often neutron densities are deter¬ 
mined with strongly interacting probes [5] such as an¬ 
tiprotons [ 311 ], proton elastic scattering |5], heavy ion 
collisions |7], Pion elastic scattering [8|, and coherent 
pion photo production [S]. Here one typically measures 
the convolution of the neutron density with an effective 
strong interaction range for the probe. Uncertainties 
in this range, from complexities of the strong interac¬ 
tions, can introduce significant systematic errors in the 
extracted neutron densities. 

It is also possible to measure neutron densities with 
electro-weak interactions, by using neutrino-nucleus co¬ 


herent scattering [IHKII] or parity violating electron scat¬ 
tering m- This is because the weak charge of a neu¬ 
tron is much larger than that of a proton. Compared to 
strongly interacting probes, parity violation provides a 
clean and model-independent way to determine the neu¬ 
tron density and likely has much smaller strong interac¬ 
tion uncertainties. In the last decades, great theoreti¬ 
cal [T2]-[T8| and experimental [M [20] efforts have been 
made to improve parity violating electron scattering ex¬ 
periments. At Jefferson laboratory, the neutron radius of 
208pb been preliminarily measured by PREX [13 na, 
and will be measured with higher accuracy by the PREX- 
H experiment |23j , while an approved experiment CREX 
aims to measure the neutron radius of "^®Ca [H]. 

In this paper, we propose to measure not only the neu¬ 
tron radius, but the full radial structure of the weak 
charge density distribution in ^®Ca, by measuring the 
parity violating asymmetry at a number of different mo¬ 
mentum transfers. This will determine the coefficients 
of a Fourier Bessel expansion of the weak charge den¬ 
sity that is model independent. By measuring the weak 
density, the full structure of neutron density can be de¬ 
rived, since the weak form factor of the neutron is largely 
known and the weak charge of the proton is very small. 
Our formalism to determine the cross-section for longitu¬ 
dinally polarized electrons scattered from "‘^Ca and the 
parity violating asymmetry Apy is presented in Section 
H. In Section HI we motivate measuring the full radial 
dependence of the weak charge density in "^®Ca and dis¬ 
cuss the large information that it contains. In Section IV 
we illustrate our formalism with an example experiment 
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and calculate the resulting statistical errors. The result¬ 
ing weak density can be directly compared to modern 
microscopic calculations of the ground state structure of 
“^^Ca using Chiral effective field theory interactions [5S] . 
We conclude in Section V that it is feasible to measure 
the full weak density distribution of ^®Ca. However this 
may be much harder for a significantly heavier nucleus 
such as ^°®Pb because more Fourier Bessel coefficients 
likely will be needed. 


II. FORMALISM 


The parity violating asymmetry for longitudinally po¬ 
larized electrons scattering from a spin zero nucleus, Apy, 
is the key observable which is very sensitive to the weak 
charge distribution. The close relationship between Apy 
and the weak charge density pw{‘>’) can be readily seen 
in the Born approximation, 

_ da/d^R +da/d^L ^ Gpq^ QwFw{<f) 
da/dflR-da/dflR 47rav^ ZFyh{q‘^) 

( 1 ) 

Here da/dQn {da/dVLi,) is the cross section for positive 
(negative) helicity electrons, Gp is the Fermi constant, 
q the momentum transfer, a the fine structure constant, 
and Fw{q^) and Fch{q^) are the weak and charge form 


factors respectively. 



Fw{q^) = j;— 
ww . 

j d^rjo{qr)pw{r) 

(2) 

Fa,{q^) = ^j 

d^rjo{qr)pyh{r). 

(3) 


These are normalized Fw(0) = F’ch(O) = 1. The charge 
density is Pch{f) and Z = j d^rpch^r) is the total charge. 
Finally, the weak charge density pwir), see Fig. and 
the total weak charge Qw = / d^rpw{r) are discussed 
below. 

The elastic cross-section in the plane wave Born ap¬ 
proximation is. 


da _ Q^cos^d) , , .2 

dn 4£;2sin'‘(f)' ^ ' ’ 


(4) 


with 0 the scattering angle. However, for a heavy nucleus. 
Coulomb-distortion effects must be included [ 12 ]. In Fig. 
I^we compare the plane-wave cross-section, Eq. to the 
cross section including Coulomb-distortion effects, see for 
example HZ). Coulomb distortions are seen to fill in the 
diffraction minima. However away from these minima 
distortion effects on the cross section are relatively small. 
The cross section calculated with the charge density from 
a relativistic mean filed model using the FSU-gold inter¬ 
action [ 25 , see Fig. agrees well with the experimental 
charge density except at the largest angles. 

We now expand the weak density of ^®Ca in a Fourier 
Bessel series. We truncate this expansion after Umax 



r(fm) 


FIG. 1: (Color online) Gronnd state (electromagnetic) 

charge and weak charge densities of '^^Ca versns radius r. The 
solid black line shows the Fourier Bessel experimental charge 
density from Ref. [I] while the dotted red line shows the 
charge density of the FSU-Gold relativistic mean field model. 
Finally the green dashed line shows the weak charge density 
of the FSU-Gold model. 


terms and assume the weak density pw{f) is zero for 
r> Rmax- This expansion will be model independent if 
truncation errors are small. 

Flmax 

Pw{r) = ^ aijo{qir) (5) 

Here qi=iTT/Rmax and jo{x) = sin(x)/x. 

To minimize measurement time we would like Umax 
and Rmax to be as small as possible while still accurately 
representing the full weak density. In this paper we con¬ 
sider 


Rmax — 7 fm, (6) 

since the weak charge density determined from many den¬ 
sity functionals is small for r > 7 fm. In addition we use, 

l^max — 6, (7) 

because the expansion coefficients determined for 
many density functional calculations of ^®Ca are very 
small for i > 6. We determine truncation errors using a 
model weak charge density based on the FSU Gold rela¬ 
tivistic mean field interaction |2S| , see below. This model 
density has 1.8 x 10“^ of the weak charge at r > Rmax = 7 
fm, and the expansion coefficients |oi| for * > Umax = 6 
are all < 7 x 10“^ fm“^. This is an order of magnitude 
or more smaller than the smallest |ai| for i < Umax- 
We now consider determining the six coefficients ai for 
i = 1 to 6. In plane wave Born approximation, a given 
Qj can be determined from a measurement of Apy(qj) at 
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FIG. 2: (Color online) Differential cross section for 2 GeV 
electrons elastically scattered from "'^Ca versus scattering an¬ 
gle. The experimental charge density is used for the solid 
black line including Coulomb distortions and the dotted blue 
line in a plane wave impulse approximation. Finally the 
dashed red line uses the model FSU Gold relativistic mean 
field theory charge density including Coulomb distortions. 


momentum transfer qj = jir/Rmax- In principle only five 
measurements are needed to determine the six because 
the weak density is normalized to the total weak charge 
/ cPrpw{r) = Qw- To be very conservative we use this 
normalization condition to determine oi. If instead we 
used the normalization to determine Og, considerably less 
beam time might be needed for a given statistical accu¬ 
racy. However, the resulting density might then be more 
sensitive to truncation errors. 

Note that in plane wave Born approximation Ap^{qj) 
is only sensitive to Oj because of the orthogonality of 
the Fourier Bessel series. When Coulomb distortions are 
included Apy{qj) is still primarily sensitive to aj and only 
depends very slightly on the other coefficients. This will 
be shown in Figs. |^and[^in Section IV. 

We consider a reference weak charge density by cal¬ 
culating pw{f) for a realistic model and then expanding 
the model density in the Fourier Bessel series. For the 
model the expansion coefficients are given by, 

_ Pw{r)3o{q^r)r^ dr 

‘ rflmax -21 "VVj” ' ■ 

Jo Jo[q^r)r^dr 

After having determined the weak density distribution, 
we can constrain the neutron density Pn(r) in ^®Ca since 
the neutron density is closely related to the weak density. 
If one neglects spin-orbit currents that are discussed in 
Ref. [26j . and other meson exchange currents m one 


can write, 

Pw{r) = J fiV'{4Gf (|r - r'|)p„(r)-b4Gp (|r - r'|)pp(r)} . 

(9) 

Here Pp{r) is the proton density and Gf (r) and G^ (r) are 
the Fourier transforms of the neutron and proton single 
nucleon weak form factors, 

4Gf (r) = Q„G|(r) + OpG^(r) - G|;(r), (10) 

4G^(r) = QpG^ (r) + Q„G^(r) - G%{t), (11) 

where G^{r) and G%{r) are Fourier transforms of the 
proton and neutron electric form factors. They are nor¬ 
malized 

J d^rGlir) = 1 , J d?TGl{r) = 0 . ( 12 ) 

Finally G’^^{r) is the Fourier transform of strange quark 
contributions to the nucleon electric form factor [55H3T] 
and is normalized / d?rG‘‘^{r) = 0. The weak form fac¬ 
tors are normalized, 

J d\ 4Gf (r) = Q„, J d\ 4G^(r) = Qp. (13) 

The weak charge of the neutron is -1 at tree level, 
while the weak charge of the proton Qp is 1 — 4 sin^ 9w 
at tree level. Including radiative corrections [32l |33] one 
has, 

Qn = -0.9878, Qp = 0.0721. (14) 

Finally, the total weak charge of “^^Ca is, 

Qw= J d^rpw{r) = iVQ„ + ZQp = -26.216 . (15) 

Further radiative corrections, for example from ^ — Z 
box diagrams [Ml |M] , are not expected to be important 
compared to this large value of Qw- 

We emphasize that parity violating experiments can 
determine the weak density pw{r) in a model indepen¬ 
dent fashion. This can be compared to theoretical pre¬ 
dictions for pw{r) that are obtained by folding theoreti¬ 
cal nucleon densities p„(r) and Pp{r) with single nucleon 
weak form factors and possibly including meson exchange 
current contributions. 


III. MOTIVATION 

In this section we discuss the information content in 
the weak charge density and some of the physics that 
would be constrained by measuring pw(r) with parity 
violating electron scattering. First the weak radius Ryj = 
[/ d^rr"^pw{r) / QwY^'^ is closely related to the neutron 












4 


radius, see for example [22]. This has been extensively 
discussed. 

The surface thickness of pwir) can differ from the 
known surface thickness of Pch{r) and is expected to be 
sensitive to poorly constrained isovector gradient terms 
in energy functionals. One way to constrain these gradi¬ 
ent terms is to perform microscopic calculations of pure 
neutron drops in artificial external potentials, using two 
and three neutron forces. Then one can fit the resulting 
energies and neutron density distributions with an energy 
functional by adjusting the isovector gradient terms. It 
may be possible to test these theoretically constrained 
isovector gradient terms by measuring the surface be¬ 
havior of pwij')- 

Next, the interior value of pw{f) for small r is closely 
related to the interior neutron density. This, when com¬ 
bined with the known charge density, will finally provide 
a direct measurement of the interior baryon density of 
a medium mass nucleus. Previously, this has only been 
extracted in model dependent ways by fitting a density 
functional to the charge density and then using the func¬ 
tional to calculate the baryon density. This is sensitive to 
the form of the symmetry energy contained in the func¬ 
tional. The interior baryon density is thought to saturate 
(stay approximately constant) with increasing mass num¬ 
ber A. This saturation density is closely related to the 
saturation density of infinite nuclear matter and insures 
that nuclear sizes scale approximately with A^/'^. 

The saturation density of infinite nuclear matter po ~ 
0.16 fm“^ is a very fundamental nuclear property that 
has proved difficult to calculate. It is very sensitive to 
three nucleon forces and calculations with only two body 
forces can saturate at more than twice po- Indeed micro¬ 
scopic calculations that use phenomenological two nu¬ 
cleon forces fit to nucleon-nucleon scattering data and 
phenomenological three nucleon forces fit to properties of 
light nuclei may not be able to make sharp predictions for 
Po because of unconstrained short range behavior of the 
three nucleon forces m- As a result these calculations 
often do not predict po but instead fit po by adjusting 
three nucleon force parameters. 

Alternatively, chiral effective field theory provides a 
framework to expand two, three, and four or more nu¬ 
cleon forces in powers of momentum over a chiral scale. 
There are now a growing number of calculations of nu¬ 
clear matter, see for example [351 [SSj- However, at this 
point it is unclear how well the chiral expansion con¬ 
vergences for symmetric matter at nuclear densities and 
above. There could still be significant uncertainties from 
higher order terms in the chiral expansion and from the 
dependence of the calculation on the assumed cutoff pa¬ 
rameter and on the form of regulators used. Furthermore, 
there are uncertainties in the energy of nuclear matter 
that arise from uncertainties in short range parameters 
that are fit to other data. Finally, some higher order 
terms in the chiral expansion may be somewhat larger 
than expected because of large contributions involving 
Delta baryons. 


In addition to calculations for infinite nuclear matter, 
advances in computational techniques have now allowed 
improved microscopic calculations directly in finite nu¬ 
clei. These calculations can directly test nuclear satu¬ 
ration by seeing how predictions compare to data as a 
function of mass number A. One can test both isoscalar 
and isovector parts of these calculations by comparison 
to both interior charge and weak charge densities. The 
interior weak density may be sensitive to three neutron 
forces and reproducing it may allow better predictions 
for very neutron rich medium mass nuclei where three 
neutron forces may also play a very important role. 

Finally we expect shell oscillations in pw('r). Shell os¬ 
cillations have been observed in Pchi^) for a variety of 
nuclei. For example there is a small increase in Pchij") for 
^°®Pb as r —>• 0 due to the filling of the 3S proton state. 
However observed shell oscillations in p^h are often much 
smaller than those predicted in many density functional 
calculations. Indeed almost all density functional calcu¬ 
lations over predict the small r bump in Pch{f) for ^°®Pb, 
see for example [ID]. It may be very useful to finally have 
direct information on shell oscillations for neutrons in ad¬ 
dition to protons. This could suggest changes in the form 
of the density functionals that are used that would cor¬ 
rect the shell oscillations. 

We conclude this section. A model independent deter¬ 
mination of pw{f) and features of the neutron density 
including surface thickens, central value, and shell oscil¬ 
lations address a number of important current problems 
in nuclear physics. Together with Pchi^) they will liter¬ 
ally provide a detailed picture of where the neutrons and 
protons are in an atomic nucleus. 


IV. SAMPLE EXPERIMENT 

In this section we evaluate the statistical error for the 
measurement of the weak charge density of ^®Ca in a 
sample experiment. As an example we consider measur¬ 
ing Apy at five points during a single run in Hall A 
at Jefferson Laboratory. The total measurement time for 
all five of the points is assumed to be 60 days. The ex¬ 
perimental parameters including beam current /, beam 
polarization P, detector solid angle AH, number of arms 
N, and the radiation loss factor C are assumed to be sim¬ 
ilar to the CREX experiment | 21 | and are listed in Table 
0 see also Ref. [H]. This example provides a conser¬ 
vative baseline for the final statistical error. Measuring 
some (or all) of the points at other laboratories such as 
Mainz, or combining data from other experiments could 
significantly reduce the statistical error. For example if 
the CREX experiment is run first and provides a very ac¬ 
curate low point, then that information could be used 
to reduce the number of future measurements needed to 
determine the full weak density. We emphasize that in 
this section we only present statistical errors. Of course 
any real experiment will also have systematic errors that 
we discuss briefly at the end of this section. 
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Parameter Value 


I 

150M 

p 

0.9 

ptar 

2.4 X 10^^ cm“2 

AU 

0.0037 Sr 

N 

2 

c 

0.34 


TABLE I: Assumed experimental parameters including beam 
current 7, beam polarization P, target thickness ptar, detector 
solid angle ATI, number of arms N, and radiation loss factor 

C- 

We calculate the statistical error in the determination 
of the Fourier Bessel coefficients of the weak density. 
The total number of electrons detected Ntot in a mea¬ 
surement time Ti is 

Ntot = mptar^C^^N . (16) 

The statistical error in the determination of is Aa^, 
— = (NtotApMifP^e^Y^^ (17) 

Oul \ / 

Here is the sensitivity of Apy to a change in and is 
defined as, 

^ _ dhiAp^[qi) _ ai dA 

pv (18) 

d In Qi Apy dtti 

In plane wave Born approximation = 1 and including 
Coulomb distortions one still has ei « 1, see Figs. |4]and 

Our reference weak charge density pwi'f’) for ^®Ca is 
the Fourier Bessel expansion of a relativistic mean field 
theory model using the FSU-Gold interaction, see Fig. 
The Fourier Bessel coefficients determined from Eq. |^are 
listed in Table II. This model yields a charge density for 
“^^Ca that agrees well with the experimental charge den¬ 
sity from Ref. [T] except in the central region as shown 
in Fig. [2 

It is crucial to use a very accurate charge density in 
the determination of the weak charge density. In Fig. 
Is] we show Apy calculated with the same reference weak 
aensity but with the FSU-Gold model or the experimen¬ 
tal charge density. There is a significant difference at 
large momentum transfers. Note both curves include 
Coulomb distortions. In the following we will always in¬ 
clude Coulomb distortions using the code ELASTC [H] 
and use the full experimental charge density. 

We now consider five measurements of Apy[qi) at mo¬ 
mentum transfers for i = 2 to 6. In general one may 
be able to improve statistics, for a given g^, by going 
to a more forward angle and higher beam energy. We 
somewhat arbitrarily restrict the scattering angle to be 
at least five degrees, since this is the scattering angle for 
the septum magnet of the PREX experiment. We also 



FIG. 3: (Color online) Parity violating asymmetry Apy for 
2 GeV electrons elastically scattered from versus scat¬ 
tering angle. The dashed red line is based on the weak and 
charge densities from the FSU-Gold relativistic mean field 
model. The solid black line uses the experimental charge den¬ 
sity from Ref. [I] and the FSU-Gold weak density. 

limit the beam energy to no more than 4 GeV as a pos¬ 
sible restriction from the HRS spectrometers in Hall A. 
Thus we consider five measurements with the kinematics 
in Table UTl 

The statistical error in the total weak charge density 
Pw{t) is the quadratic combination of the errors for the 
six Fourier Bessel terms. Note that the statistical errors 
Aoi for different i are independent. 

Apu,eafc(r) = I ^ |^Aaijo(g*r)j j"" (19) 

i=l 

The individual errors Aa^ depend on the time Ti spent 
measuring Apy{qi) at momentum transfer g^. As a sim¬ 
ple example we optimize the individual Ti, subject to a 
constraint on the total measurement time 
6 

^ T, = 60 days, (20) 

i=2 

in order to minimize the statistical error in pw{A) at 
r = 0. Note that the error in oi is calculated by using 
the normalization condition f pwr(r)(Pr=Qw■ The indi¬ 
vidual Ti and the fractional errors in each of the Aa^ are 
listed in Table El 

Table[lT|shows that most of the time is spent measuring 
the highest momentum transfer points. This is because 
the cross section falls so rapidly with increasing g. One 
alternative, to this model independent approach, would 
be to constrain the higher i coefficients ai from theory 
and only measure Apy for smaller momentum transfers. 
This could significantly reduce the run time and the sta¬ 
tistical error. 
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Qi 

fin"^ 

E 

GeV 

dcT 

dQ. 

mb 

■^pv 

ppm 

T 

days 

ai 

fm”^ 

Aaijui 

% 

0.45 





0.0752 

1.1 

0.90 

2.06 

2.44 

2.54 

5 

0.0468 

5.9 

1.35 

3.09 

1.07 X 10~^ 

8.31 

7 

-0.0438 

7.6 

1.80 

4 

2.9 X 10"® 

9.92 

10 

-0.0147 

27 

2.24 

4 

4.05 X 10“^ 

22.5 

15 

0.0161 

29 

2.69 

4 

9.7 X 10"® 

36.5 

23 

0.0066 

90 


TABLE II: The momentum transfer qi, beam energy E, cross 
section, parity violating asymmetry measnrement time 
T, Fourier Bessel expansion coefficient ai of the weak charge 
density as determined for the FSU-goId relativistic mean field 
model, and fractional statistical error Aoi/fli. Note that the 
error in ai is determined by normalizing the weak charge den¬ 
sity to the total weak charge Qw- 



FIG. 4: (Color online) Parity violating asymmetry Ap^ for 
2 GeV electrons elastically scattered from ^®Ca versus scat¬ 
tering angle. The solid black line shows the asymmetry curve 
based on the charge density from Ref. [I] and the FSU-Gold 
weak density. The red dotted, green dashed, and blue dot- 
dashed curves show Ap^ when 02 , 0,3 and 04 have been varied 
one at a time by ten percent respectively. The arrows show 
the momentum transfers qi for i = 2, 3 and 4. 


Figure shows the statistical error in the weak density 
Pw{f) as a function of radius r. The error in pw is largest 
for small r and gradually decreases as r increases. Thus 
it is most difficult to determine pwi^) near the origin. 
There may be several ways to decrease the error band in 
Fig. § One could measure with higher beam currents 
and or with larger acceptance spectrometers. Alterna¬ 
tively, one could measure for a larger time either as one 
extended experiment or by combining experiments that 
each focus on only some of the qi points. 

We now discuss systematic errors. Because it is so dif- 
hcult to get good statistics for large q, the higher i coef- 
hcients ai may only be determined with somewhat large 



FIG. 5: (Golor online) As per Fig. I^but for larger scattering 
angles. The arrows show qi for i = 4, 5 and 6. 



FIG. 6: (Golor Online) Weak charge density Pwij) of ‘*®Ca 
versus radius r. The black solid line is the reference FSU-Gold 
weak density and the black dotted lines show the statistical 
error band that could be obtained by measuring at hve 
specific momentum transfers with a total running time of 60 
days, see text. The dashed red line shows the experimental 
(electromagnetic) charge density [1]. 


statistical errors Aa^/ai. As a result many systematic 
errors such as determining the absolute beam polariza¬ 
tion or from helicity correlated beam properties may be 
less important. Instead backgrounds, from for example 
electrons that scatter from collimators used to define the 
acceptance, could be important because the elastic cross 
section is small (at higher q). 

We have focused on ^®Ca. Determining the full pw{f) 
for a significantly heavier nucleus such as ^°®Pb may be 
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dramatically harder. This is because more Fourier Bessel 
coefficients will likely be needed and because the cross 
section drops extremly rapidly with increasing q. Thus 
it may be very hard to measure Apy, at high enough 5 , 
in order to directly determine the weak density in the 
center of ^°®Pb. 


V. CONCLUSIONS 

The ground state neutron density of a medium mass 
nucleus contains fundamental nuclear structure informa¬ 
tion and it is at present relatively poorly known. In this 
paper we explored if parity violating elastic electron scat¬ 
tering can determine not just the neutron radius, but the 
entire radial form of the neutron density Pn{f) or weak 
charge density pw{'>') in a model independent way. We 
expanded the weak charge density pw(.f) in a model in¬ 
dependent Fourier Bessel series. For the medium mass 
neutron rich nucleus “^^Ca, we find that a practical parity 
violating experiment could determine about six Fourier 
Bessel coefficients ai and thus deduce the full radial struc¬ 
ture of both pw(.f) and the neutron density p„(r). The 
resulting pw (r) will contain fundamental information on 


the size, surface thickness, shell oscillations, and satura¬ 
tion density of the neutron distribution. 

Future work could optimize our model experiment to 
further reduce the statistical errors by for example using 
large acceptance detectors and combining information 
from multiple experiments and or laboratories. Future 
theoretical work exploring the range of weak charge den¬ 
sities to be expected with reasonable models and micro¬ 
scopic calculations would also be very useful. The mea¬ 
sured pw{t)^ combined with the previously known charge 
density Pchix), will literally provide a detailed textbook 
picture of where the neutrons and protons are located in 
an atomic nucleus. 
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